AD/A-004  827 


THEORY  AN D C A LC U LAT  ION  FOR  DISPLACE- 
MENTS AND  STRESSES  AT  AN  EARTHQUAKE 
SOURCE 

Paul  G.  Richards 

Lar.iont  - Done  rty  Geological  Observatory 


Pre  pare  d for : 

Air  Force  Office  of  Scientific  Re  search 
Advanced  Res  6;  arch  Projects  Agency 


30  October  1974 


DISTRIBUTED  BY: 


National  Technical  Information  Service 
U.  S.  DEPARTMENT  OF  COMMERCE 


UNCLASSIFIED 


SECURITY  CL  ASSIFu.  i ' OF  THIS  P»OC  'i  V.n  P.!S  Fnlueifl 


Kl'.AI)  MSTRUC 


REPORT  DOCUMENTATION  PAGE 


BKFOKK  COMPIJC'i 


3.  RLCIPILMT'S  CATALOG 


2.  GOV1  ACCESSION  KO 


(&>/  AjlWA/a 


S.  TYf't  or  REPORT  A PERIOD  COVi 

Semi-Annual,  Scienti 
interim 


4.  TITLE  (end  Subtitle) 


THEORY  AND  CALCULATION  FOR  DISPLACEMENTS 


6.  performing  org.  report  number 


and  STRESSES  AT  AN  EARTHQUAKE  SOURCE 


contract  or  grant  number^*. 


7.  author^; 


pi]  11620-7*4-0-0029 


Professor  Paul  G.  Richards 


10.  PROGRAM  U r.MEllT.rr-.  >:  -T.  TASa 


PERFOniVKO  ORGANIZATION  name  and  address 

Columbia  University 

Lainont-Doherty  Geological  Observatory 
Palisades,  Mew  York  1096-1 


ARIA  & MiOKk  UNIT  NUM. 

ARPA  Order  Mo. 
Program  Code  Mo 
Program  Element 


12.  REPORT  DAT  E 

30  October 


CON  1 i.  LLINC  OFFICE  NAME  AND  AODRfcSS 

Advanced  Research  Projects  Agency 
1401  Wil.son  Eoulevard 

Arlington,  Virrjnia  ?.??.  

MONITORING  AGENCY  NAME  e AOORESf  il  different  Row  Controlling  Office) 


13.  NUMBER  OF  PAGES 


Air  Force  Office  of  Scientific  Research 
1400  Wilson  Boulevard  '' 
Arlington,  Virginia  22209 


10.  DISTRIBUTION  ST  ATI  MFMT  (o!  thl-  ' .r  "> 


Ar\rTr»oved.  for  nubile  release 


*7.  DISTRIBUTION  STATEMENT  (ol  the  abstrect  entered  In  Hod.  20,  If  different  Ira-  l<eport) 


16.  SUPPL  EMENT  ARY  NOT  ES 


Reprodwcod  by 

NATIONAL  TECHNICAL 
INFORMATION  SERVICE 

US  D*p»r«men»  ol  Commerce 
Spr.ngf.0ld,  VA.  2215* 


19.  KEY  WORDS  (Continue  on  reverse  tide  if  nece  t sery  nnd  Identity  by  block  number) 


Earthquake  - explosion  discrimination 


ourco  theory 


Earthquake 


20.  ADST  PACT  (Continue  on  reverse  ride  If  neces  secy  end  Identify  ty  block  number) 

This  progress  report  describes  the  near-field  motions  of  an  oart 
quake  source  model  which  is  governed  by  simple  physical  constrai 
in  terns  of  dynamic  friction  on  a fault  plane.  It  also  describes 

comnuti  rorrai  which  will  sulfite  : ’ . ril : >1  *??  :j 


EDITION  OF  I IIOV  t‘.  IS  ODSOLri  E 


UMCl.AKSIFJKD 


H*>»1 
O W 


ARPA  Order  Number: 

Program  Code  Number: 
Contractor: 

Effective  Date  of  Contract: 
Contract  Expiration  Date: 
Amount  of  Contract: 

Contract  Number: 

Principal  Investigator: 
Program  Manager: 

Title  of  Work: 


1827-11 

^FIO 

Columbia  University 
1 October  1973 
30  September  1975 
$70,000.00 
F M620-74-C-0029 
Paul  G.  Richards,  91^-359-2900 
Paul  G.  Richards,  91ii-359-2900 
Displacements  and  Stresses  for 
an  Earthquake  Source 


1. 


SUMMARY 


During  the  second  six  months  of  the  subject  contract,  the 
purpose  of  the  research  program  was  to  investigate  properties  of 
an  earthquake  source,  to  provide  a physical  basis  for  discrimination 
between  different  types  of  seismic  events.  The  two  main  approaches 

used  were : 

1.  To  investigate  thoroughly  the  only  known  analytic  solution 
for  a three  dimensional  problem  of  earthquake  faulting,  in  which 
shear  stresses  on  the  moving  parts  of  the  fault  surface  are  deter- 


mined by  a simple  law  of  friction. 

2.  To  develop  the  computational  capability  co  determine  the 

seismograms  which  would  result  from  any  assumption  about 

the  relative  shearing  motions  taking  place  across  an  earthquake 

fault . 

Research  on  the  first  of  these  approaches  is  based  on  modeling 
the  rupture  process  of  earthquake  faulting  by  a sh_ar  crack, 
shear  stress,  on  moving  faces  of  the  crack,  is  taken  as  a constant 
(the  dynamic  coefficient  of  friction)  times  the  normal  stress 
acting  across  the  crack.  The  crack  is  presumed  to  grow  from  a 
point  (i.e.,  the  point  at  which  rupture  begins  on  an  earthquake 
fault),  and  then  grows  steadily  as  an  ellipse  with  fixed  eccentricity. 
The  great  virtue  of  this  earthquake  model  is  that  it  is  based  on 
very  simple  physical  principles  of  mechanical  failure,  and  no 
a priori  assumption  is  made  about  the  displacement  on  the 
surface  Itself.  The  resulting  accelerations  can  be  found  everywhere 
within  the  elastic  medium.  Principal  properties  of  the  computed 
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solution  are:  (i)  that  particle  velocities,  near  the  rupture 

front,  increase  with  increasing  distance  from  the  point  of  initial 
rupture;  (ii)  that  it  is  exceedingly  difficult  to  detect  any  fea- 
tures in  a record  section  of  near-field  seismograms  which  v/ill 
enable  one  to  determine  the  rupture  velocity;  (iii)  that  the  con- 
cept of  "maximum  acceleration",  observed  in  t)-  e near  field  for  a 
given  earthquake,  is  an  unstable  property  of  the  positioning  of 
the  observing  point  within  the  radiation  pattern,  rather  than  a 
macroscopic  property  of  the  earthquake  itself. 

The  goal  of  our  second  major  research  approach  has  been  fully 
attained.  That  is,  we  have  developed  a computer  program  to  find 
the  seismograms  .vhich  would  result  from  virtually  any  earthquake 
source,  described  in  terns  of  the  rupture  velocities  and  relative 
displacements  taking  place  on  a fault  plane  within  an  infinite 
homogeneous  elastic  medium.  We  have  used  this  program  extensively 
to  examine  the  ground  motion  which  would  take  place  near  an  earth- 
quake fault  for  several  dislocation  models  of  the  earthquake. 

Quite  different  dislocation  models  can  give  very  similar  displace- 
ments at  stations  located  only  one  fault  dimension  from  the  fault. 
In  particular:  1)  models  with  fault  motions  and  rupture  geometry 

quite  different  from  Haskell's  propagating  ramp  can  often  have 
near-field  motions  very  similar  to  those  from  a ramp  model;  and 
2)  quite  similar  near-field  motions  are  found,  for  different  ramp 
models  in  which  the  rupture  velocity  and  rise  time  are  varied 
together  over  rather  large  ranges  (e.g.,  a factor  of  three  in 
rise  time ) . 


We  infer  that  there  is  considerable  ambiguity  in  interpreting 
near-field  displacement  data  in  terms  of  a model  of  fault  motion. 

The  ambiguity  may  be  reduced  by  using  acceleration  data,  although 
this  entails  a considerable  increase  in  computer  time  for  solving 
forward  problems. 

Both  of  our  research  approaches  indicate  a result  which  will 
be  of  importance  in  determining,  from  seismic  data  taken  in  the 
very  near-field  of  a propagating  fault,  details  of  the  fault  motion 
itself.  We  have  found  that  the  component  of  (vector)  displacement 
which  is  most  sensitive  to  the  fault  motion  is  the  component  per- 
pendicular to  the  fault  surface.  One  might  have  expected  that 
the  most  sensitive  component  would  be  that  parallel  to  the  direction 
of  slip  of  the  fault.  Many  experiments,  particularly  rock  mecha- 
nics experiments,  have  emphasized  measurement  of  fault  motion 
parallel  to  the  slip  direction,  but  we  find  this 

is  a component  which  is  very  insensitive  to  rupture  velocity. 

The  tv.'o  research  approaches  described  above  are  each  expected 
to  lead  to  a research  publication.  Our  efforts  in  the  tnird  six- 
month  period  of  the  contract  will  be  directed  tov.’ards  problems 
in  which  the  free  surface  of  the  earth  has  an  important  effect 
on  waves  radiating  from  an  earthquake  source. 

Elastic  Waves  near  the  Rupture  Front  of  a Growing  Plane  Elliptical 

Shear  Crack: 

A Three-Dimensional  Solution 

Earthquakes  are  a stress  relaxation  phenomenon:  upon  this, 

there  has  been  general  agreement  among  seismologists 


r over  fifty 


u. 

years.  In  this  section,  we  describe  a kind  of  fault  motion  which 
falls  into  this  category,  in  which  an  initially  high  stress  field 
is  relieved  by  motion  on  the  fault  surface.  We  then  describe  the 
accelerations  and  displacements  in  the  near-field  of  such  a propa- 
gating fault. 

We  suppose  that  an  infinite  elastic  medium  is  subjected  to 
tectonic  shearing  stresses  in  the  xz  direction,  with  a fault  in 
the  plane  z = 0.  In  this  situation,  a crack  nucleates  at  a point, 
and  a rupture  front  spreads  steadily,  initiating  relative  motion 
across  the  fault  plane.  We  suppose  that  the  rupture  front  main- 
tains the  shape  of  an  ellipse,  shown  in  Figure  la,  spreading 
steadily  with  speed  a in  the  x direction  and  speed  v in  tho  y 
direction.  This  much  of  the  description  of  earthquake  faulting 
is  "kinematic".  However,  the  boundary  condition  we  impose  on 
fault  plane  is  "dynamic".  We  suppose  that  dynamic  friction 
forces  act  on  those  parts  of  the  fault  which  are  in  relative 
motion:  shearing  stresses  are  taken  to  be  proportional  to  the 

normal  stress.  It  is  often  convenient  to  use  spherical  polar 
coordinates,  shown  in  Figure  lb,  rather  than  cartesian  co- 
ordinates. Our  computations  are  principally  for  points  near 
the  fault  plane,  since  we  are  interested  in  trying  to  detect  t 
the  presence  of  the  rupture  front  as  it  moves  past  observing 
stations.  Many  of  our  computations  are  carried  out  for  the 
case  0 = 85°.  The  rupture  speed  is  taken  as  less  even  than  the  shear 
wave  speea,  so  the  only  wave  fronts  in  this  problem  are  tho 
spherical  P and  S wave  fronts  emanating  from  the  point  of  initial 
nucleation.  The  rupture  front  itself  does  not  radiate  a wave 
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front . 

These  assumptions  turn  out  to  be  enough  to  set  up  a boundary 
value  problem  for  motions  throughout  the  elastic  medium.  The 
solution  has  been  published  by  Richards  (1*73),  and  here  w?  are 
concerned  with  evaluating  the  solution  to  obtain  insight  into  the 
behavior  of  shear  stresses  and  accelerations  in  the  vicinity  of 
a propagating  fault. 


Figure  2 shows  one  of  the  quantities  which  can  easily  be 
computed.  On  the  bottom  left  is  shown  a quadrant  of  the  elliptical 
part  of  the  fault,  which  is  in  motion.  In  this  case,  rupture  in 
the  x direction  is  taken  as  the  Rayleigh  speed,  ar-»  in  the  y 
direction  as  the  shear  wave  speed..  The  xz  component  of  stress 
is  being  relieved  by  the  fault,  and  the  major  part  of  the  figure 
shows  four  different  time  histories  for  the  stress  rate  r „ , at 
different  azimuthal  angles  4>,  for  points  near  the  fault  plane. 

It  is  the  stress  rate  which  comes  naturally  out  of  the  solution 
method,  and  it  is  natural  to  express  the  solution  in  a dimension- 
less form.  The  factors  relating  dimensionless  and  dimensioned 
stress  rates  and  times  are  shown  at  the  upper  right  of  Figure  2: 
a,  B and  p are  respectively  the  P wave  speed,  the  shear  wave  speed, 
and  the  density;  a is  the  particle  velocity  at  the  point  of 
nucleation.  Particular  values  of  these  dimensioning  factors 
are  shown,  for  a point  10  km  away  from  the  initial  rupture, 
with  an  assumed  particle  velocity  of  10  cm  per  second  the 

initial  rupture. 
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Figure  3 shows  the  result  of  time  integrating  Figure  2,  to 
obtain  rX2  itself.  Again,  the  shear  stress  is  plotted  in  dimen- 
sionless form  for  four  points  near  the  fauxt  plane,  differing 
only  in  azimuthal  angle . It  may  be  seen  that  txz  has  a rather 
different  time  history  at  each  of  the  four  azimuths,  but  the 
value  for  large  times  is  everywhere  the  same,  about  -8  units. 

With  an  initial  particle  velocity  of  100  cm  per  second,  this 
translates  to  a stress  drop  of  about  200  bars.  Figure  <1  differs 
from  figure  3 only  in  that  the  rupture  velocities  are  halved, 
lor  the  same  initial  particle  velocity  of  100  cm  per  second, 
the  stress  drop  now  is  about  300  bars. 

We  should  like  to  emphasize  the  three-dimensional  nature 
of  this  earthquake  model.  Because  the  xz  component  of  stress 
is  the  only  component  being  relieved  by  fault  motion,  it 
follows  that  the  other  shear  component  on  the  fault  plane, 

Tyz,  should  have  final  values  unchanged  from  zero.  Figure  5 
shows  this  indeed  to  be  the  case,  in  that  xyz  at  different  posi- 
tions near  the  fault  plane  does  tend  to  zero  for  large  times, 
but  it  will  be  seen  at  Intermediate  values,  particularly  at  times 
near  the  rupture  arrival,  to  be  large.  This  kind  of  feature 
is  important  to  consider,  since  fault  motions  must  begin  at  a point, 
and  then  spread  over  the  fault  surface  with  a rupture  front  which 
is  a closed  curve  (or  several  closed  curves)  in  the  fault  plane. 

Figures  2,  3,  4 and  5 all  describe  the  stress  history  near 
a rupturing  fault  and  help  us  understand  how  stress  is  relieved 
by  the  stress  relaxation  phenomenon  of  the  earthquake.  But,  of 
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course,  seismic  instruments  are  built  typically  to  record  accelera- 
tion, or  displacement c Figure  6 shows  both,  again  for  four  differ- 
ent azimuths,  and  plotted  here  as  solid  lines  are  the  x components 
of  acceleration  and  displacement,  for  points  near  the  fault  plane. 
The  dashed  line  is  the  displacement  time  history  for  that  point 
on  the  fault  surface  which  is  nearest  the  station  with  azimuthal 
angle  as  shown,  and  8 = 85°.  Acceleration  clearly  has  a different 
time  history  at  these  locations,and  if  particle  velocity  is  100  cm 
per  second  on  the  fault,  then  maximum  accelerations  are  around 
1 g at  the  distance  10  km,  for  $ = 1°,  30°,  and  60°.  The  case 
<t>  = 89°  is  pathological,  because  there  the  rupture  is  essentially 
riding  on  the  S wave  wavefront,  and  accelerations  are  up  around 
several  tens  of  g,  due  to  the  Doppler  focusing  of  energy.  However, 
the  displacement  (obtained  by  double  integration  of  acceleration) 
is  about  the  same  at  all  four  different  positions, and  in  each  case 
the  computed  displacement  is  similar  to  the  dashed  curve,  the 
actual  shear  dislocation  of  the  nearest  point  on  the  fault  surface. 

Another  displacement  that  is  important  to  investigate  is  the 
component  normal  to  the  fault  plane.  Figure  7 shows  the  accelera- 
tion of  displacement  in  this  direction,  in  the  case  of  somewhat 
slow  rupture  propagation.  The  first  remarkable  point  is  that  the 
acceleration  time  history  has  essentially  the  same  shape  at  all 
the  four  stations.  However,  the  scale  is  very  different.  This 
similarity  in  shape  is  due  to  the  fact  that  that  part  of  the  fault 
plane  which  is  in  motion  has  been  twisted  through  a constant  angle 
out  of  its  original  plane,  so  that  the  z components  of  acceleration 
and  displacement  will  have  a magnitude  which  is  proportional  to  the 
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x coordinate.  Displacement  is  shown  here  as  a dashed  line.  The 
second  important  point  to  make  is  that  the  z displacement  may  be 
a good  marker  of  the  rupture  front  going  nearby.  If  one  actually 
had  data  for  displacement  of  the  type  shown  as  dashed  lines  in 
Figure  7,  then  one  would  not  be  far  wrong  in  picking  the  times 
of  rapid  change  in  z displacement  as  the  time  the  rupture  front 
is  passing  nearby. 

The  stress  relaxation  model  being  described  in  this  section 
has  the  property  of  a stress  drop  which  is  proportional  to  particle 
velocity  at  the  initial  point  of  rupture.  It  is  clearly  useful 
to  know  the  constant  of  proportionality,  and  in  Figure  8a  is 
plotted  a value  for  this  constant,  for  a range  of  rupture  velocities. 
The  solid  curve  in  Figure  8a  was  attained  by  evaluating  an  integral 
given  by  Kostrov  (196*0,  for  a circular  crack . Another  way  to  get 
this  curve  is  to  take  the  time  history  of  stress  rate  as  shown 
in  Figure  2 for  a point  near  the  fault  plane,  and  to  integrate 
this  up  to  find  the  final  static  value  of  the  stress  drop.  Carrying 
this  procedure  out  at  two  different  rupture  speeds  results  in  the 
two  crosses  shown  on  Figure  8a:  the  agreement  with  Kostrov’s 

curve  gives  one  confidence  that  the  numerical  work  is  self  con- 
sistent. The  basic  observation  to  make  about  Figure  8a  is  that 
the  curve  is  almost  a straight  line  through  the  origin.  In  fact, 
the  slight  negative  curvature  is  important  , since  it  implies  that 
faults  are  unlikely  significantly  to  overshoot  their  final  static 
displacement.  One  can  conclude  this  by  re-plotting  the  curve, 
using  it  to  find  out  how  much  displacement  there  is  at  the  center 
of  the  fault,  at  the  time  the  fault  surface  has  acquired  a fixed 
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radius,  e.g.,  1 km.  If  rupture  were  then  to  stop,  we  ask:  "Would 
displacement  across  the  fault  continue  to  grow  in  the  same  direction, 
or  would  it  have  to  reverse?"  The  curve  shown  in  Figure  8b  in- 
dicates that  if  the  fault  displacements  are  to  end  up  in  a final 
static  configuration  which  is  appropriate  to  the  stress  drop  pre- 
vailing during  dynamic  rupture,  then  the  displacement  will  have  to 
grow  up  to  the  static  value  from  below,  and  overshooting  is  not 
required . 
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Figure  2 


Figure  3 

Figure  4 
Figure  5 
Figure  6 


FIGURE  CAPTIONS 


Parameters  for  a growing  elliptical  crack:  (a)  the 

plane  z = 0,  seen  from  the  side  with  z > 0;  (b)  the 
definition  of  spherical  polar  coordinates  (R,  0,  (j>). 
Azimuth  angle  <f>  is  taken  between  the  plane  y = 0, 
and  the  plane  containing  the  field  point  and  the 
z-axis . 

The  time  history  of  stress  rate,  shown  at  four 
different  azimuths.  At  bottom  left  is  shown  the 
relative  positions  of  that  part  of  the  fault  surface 
which  is  in  motion  (shown  shaded);  the  S wave 
wavefront;  and  the  P wave  wavefront.  Stress  rate 
and  time  are  plotted  as  dimensionless  quantities, 
and  factors  shown  at  the  upper  right  may  be  used 
to  convert  the  plots  to  dimensioned  values  for 
specific  cases. 

As  for  Figure  H but  showing  stress.  This  figure 
was  obtained  via  a time  integration  of  Figure  2. 

As  Figure  3,  but  with  rupture  speeds  halved. 

As  Figure  3,  but  for  the  TyZ  component  of  shear. 

As  Figure  2,  but  showing  the  x-component  of  acce- 
leration and  displacement.  Acceleration  values  are 
obtained  from  the  peaked  curves,  and  left-hand 
vertical  scales:  displacements  from  the  monotonic 

curves,  and  rxght-hand  scales.  Dashed  curves  show 
the  displacement  time-history  at  nearest  points  on 
the  fault  surface  itself. 
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Figure  7 


Figure  8 


As  Figure  6,  but  for  the  slower  rupture  speeds,  and 
the  z-component  of  acceleration  (solid)  and  dis- 
placement (dashed). 

(a)  Evaluation  of  formulae  given  by  Kostrov  (196*0 
for  the  particle  velocity  at  fixed  stress-drop, 
shown  as  a function  of  rupture  speed,  for  a circular 
shear  crack. 

(b)  Fault  displacement  (i.e.,  half  the  maximum  offset), 

shown  as  a function  of  rupture  speed,  at  the  time 

* 

the  fault  acquires  a radius  of  1 kilometer. 
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Program  to  Compute  Seismograms  from  Prescribed  Shear  Dislocations : 
Comparison  of  Strong  Ground  Motion  from  Several  Dislocation  Models 

• * 1 i i ■■  ■ - . — ■ — 

One  approach  to  studying  details  of  the  earthquake  source 
mechanism  is  to  study  records  of  strong  ground  motion  near' 
the  causative  fault.  At  these  short  distances,  where  the 
effects  of  scattering  and  attenuation  are  minimal,  one  can 
best  expect  to  learn  about  the  way  rupture  grows  and  stops, 
and  the  time  function  of  displacement  on  the  fault.  However, 
interpretation  of  strong  motion  accelerograms  recorded  near 
the  causative  fault  is  not  straightforward.  Several  workers 
have  studied  these  records  in  terms  of  a dislocation  model 
in  an  infinite  space.  For  the  time  history  of  dislocation, 

Aki  (1963)  used  a propagating  step,  but  later  workers  (for 
example,  Haskell,  1969;  Kanamori , 1972;  Yrifunac,  197^; 

Trifur.ac  and  Udwadia,  1 9 7 ; Anderson,  1 9 7 zO  used  a propagating 
ramp.  These  workers  found  particular  combinations  of  para- 
meters such  that  the  displacements  resulting  from  their  models 
approximately  fit  the  integrated  accelerograms  for  the  earth- 
quakes studied.  Their  results  demonstrate  that  the  long  period 
components  of  near  field  data  can  be  fairly  well  matched,  in 

spite  of  the  limitations  of  the  theory  and  the  dislocation 
model . 

t 

Dislocation  models  are  generally  found  by  trial  and  error. 
The  model  found  is  not  necessarily  the  best,  nor  is  it  the  only 
one  which  will  adequately  fit  the  data.  In  fact,  models  signi- 
ficantly different  from  each  other  are  possible.  This  ambiguit; 
is  demonstrated  by  the  variety  of  solutions  from  several  studio: 


70 

of  the  Parkfield,  California,  earthquake  of  June  27,  1966. 
For  example,  models  have  been  found  with  rupture  velocities 
of  2.2  km/sec  (Aki,  1968;  Tsai  and  Patton,  1972),  2.4  to  2.5 
km/sec  (Trifunac  and  Udwadia,  1974)  and  2.0  to  3.0  km/sec 
(Anderson,  1974). 

In  this  paper  the  nature  of  this  ambiguity  is  examinee 
in  more  detail.  This  has  been  done  in  the  follovjing  ways. 

1)  A computer  program  has  teen  developed  which  can 
calculate,  for  a wide  variety  of  dislocation  mo- 
dels, the  time  history  of  displacement  at  any 
point  in  an  infinite,  homogeneous,  elastic  space. 

2)  Several  type.-,  of  dislocation  moaels,  including  a 
propagating  ramp,  are  used,  and  displacements  at 
nearby  stations  are  studied  to  try  to  learn: 

A)  What  factors  of  the  fault  motion  are  important 
to  ground  displacement  nearby,  and  B)  From  the  near 
ground  displacement,  how  much  can  be  learned  about 
details  of  the  fault  dislocation. 

3)  Propagating  ramp  models  permit  relatively  rapid 
computation,  so  we  have  used  them  in  several  cases 
to  compare  their  near-field  motions  with  those  of 
otner  models.  A surprisingly  good  match  can  often 
be  found . 

4)  An  attempt  is  made  to  understand  why  a propagating 


I 


ramp  can  often  work  so  well  in  fitting  displacements 
caused  by  other  dislocation  models. 


*2/ 


Description  of  Method 

We  model  the  earthquake  by  a discontinuity  in  displace- 
ment across  a rupturinc  fault.  Specifically: 

1.  A. fault  plane  is  assumed,  here  the  x3  = 0 plane. 

A location  on  the  fault  plane  is  given  by 

£ = C2K  Station  locations  (observation 

points)  are  given  as  x = (x^  x2 , x3). 

2.  Here  the  fault  will  refer  to  that  bounded  region 

of  the  fault  plane  which  is  offset  during  the 

* 

earthquake.  During  the  earthquake,  as  viewed 

in  inertial  space,  a particle  at  point  C attached 
to  the  x3  > 0 side  of  fault  has  the  motion  u (£,  t). 

If  the  same  particle  were  attached  to  the  x3  < 0 

side  of  fault,  it  would  have  the  displacement 

u (£,  t).  The  dislocation  function  is  then 
Ut,  t)  = uUl,  t)  - U*(t,  t). 

3.  Using  formulae  given  by  Haskell  (1969),  the  dis- 
placement at  any  point  in  an  infinite  homogeneous 
elastic  space  caused  by  the  dislocation  D(£,  t) 
can  be  calculated. 

A computer  program  was  developed  which  is  capable  of 

4-  4 

using  any  dislocation  function  D(£,  t)  provided  that  function 
Is  smooth  enough  in  space  to  be  integrated  using  the  ROMBERG 

scheme.  This  differs  from  the  application  usually  made  of 

Haskell’s  formula,  in  that  we  do  not  assume  D to  be  the  pro- 
duct of  a function  of  space  and  a function  of  time.  Furthermore, 

we  do  not  restrict  the  time  dependence  to  be  ramplike.  V/e  do 


however,  require  that  the  function  be  continuous  in  tine. 
A typical  formula  from  the  Haskell  paper  is 


u 


t ) 


( $2/^tt  ) 


Ca, 


t-t 1 ) t 1 dt 1 


+2(6yj  2-l) (ar)“2  D(S,,  $2, t-r/a) 


-3(^Y,2-l)(&r)-2  DU,,  £2  ,t-r/3) 


+2Y12(a3r)"1  DU,,  52,  t-r/a) 


~(2Yi2-l)(C3r)“l  DUi,  Sa,  t-r/G) J d£ ,dC2 
Here  ui  is  the  displacement  in  the  x,  direction  due  to  a dis- 
location 5 = (D,  0,  0)  on  the  fault  plane  Z.  Parameters  a,  8 
are  the  P-  and  d-  wave  velocities;  r = |x  - £| ; and 
Yi  * (*i  - £i)/r. 

D 

Integration  of  the  terms  in  r were  achieved  by  reversing 


the  order  of  integration  and  differentiation  as  follows: 
YjY,  2 ( a 3 r ) “ 1 DUi,  Ca,  t-r/a' 


“ Y3  (2Yi  2-l)  (33r)-1  DU,,  C2,  t-r/{J) 


dS: 


= dt  l|l2Y3Y,2(a3r)-1  DU,,  C2,  t-r/a) 


( 2 Y , 2-l) (B'r)”1  D(5j , C2,  t-r/G)j d£ ,d£2 


For  this  step  to  be  valid,  it  is  usually  sufficient  that  the 
dislocation  function  be  continuous  in  time.  This  allov/s  the 
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use  of  dislocation  functions  which  have  infinite  particle 
velocities  at  the  rupture  front.  The  function  which  results 
from  the  double  integration  over  the  fault  plane  varies  smoothly 
in  time  and  does  not  cause  any  unusual)  problems  in  numerical 
differentiation . 

All  the  ntegration  is  done  using  a ROMBERG  scheme.  The 
ROMBERG  integrations  were  stopped  when  two  successive  esti- 
mates of  the  integral  differed  by  less  than  a specified  er-or, 
usually  taken  as  5%.  At  worst,  the  program  took  about  V*. 
sec/Jata  point  (3  displacement  components  for  the  sTame  point 
in  time)  on  an  IBM  360-95,  and  at  best  was  ten  times  faster. 

The  program  was  chocked  in  the  following  ways: 

1.  Using  a propagating  ramp  dislocation  function,  dis- 
placement seismograms  were  found  to  agree  well  with 
results  of  Tsai,  Trifunac,  aid  Boore  (personal 
communications ) . 

2.  For  several  models,  calculation  of  displacements  at 
stations  located  on  planes  of  symmetry  resulted  in 
displacement  amplitudes  (for  theoretically  zero  com- 
ponents) at  least  five  orders  of  magnitude  si  aller  than 
for  non-zero  components. 

3.  Using  a propagating  ramp  dislocation  function,  an 
array  of  19  points  in  space  for  three  times  was  found 
to  yield  a finite  difference  form  of  the  elasticity 
equation : 

£ ' 
u « (a2  - + fi2V2u 
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A severe  test  was  made  by  taking  the  dislocation  D 
to  be  that  for  a self  similar  crack,  and  comparing 
the  calculated  displacements  with  results  obtained 
using  an  entirely  different  method  (Richards,  1973). 

The  dislocation  functions  used  here,  which  are  intended 
to  be  illustrative  and  are  not  physical]y  realistic,  are  given 

. i 

belovi . In  Fig.  1,  an  example  of  all  but  one  of  the  dislocation 
functions  is  given  to  show  how  it  grows  in  space.  The  shading 
(see  caption)  denotes  either  that  part  of  the  fault  which  is 
in  motion,  cr  .hat  part  which  has  achieved  less  than  half  tne 
final  average  offset.  The  symbol  given  for  each  dislocation 
model  in  Fig.  1,  which  schematically  suggests  the  way  the  faults 

grow,  will  be  used  to  identify  models  in  later  figures. 

The  self-similar  crack  is  a solution  (Burridge  and  Willis, 

1969)  for  a crack  growing  in  an  infinite  anisotropic  elastic 
medium  with  constant  tractions  on  the  fault  surface.  To  be 
useful  for  our  purposes , there  must  also  be  a mechanism  for  stop- 
ping the  rupture.  Thus  tie  dislocation  function  here  called 
v "seif-similar  crack"  Is,  initially,  a self' similar  crack, 

but  it  is- stopped  instantaneously.  The  dislocation  function 

✓ 

called  a "pinned  annealing  crack"  also  begins  as  a self-similar 

crack.  When  the  rupture  front  reaches  a designated  point  on 
* 

the  fault  plane,  that  point  does  not  rupture.  From  that  point, 
an  annealing  iront  spreads  which  stops  the  rupture  instanta- 
neously as  it  passes  other  points  on  the  fault. 


1 


^5~ 


The  one  source  omitted  from  Fig.  1 is  a smoothed  propa- 
gating ramp  with  the  time  function  suggested  by  Ohnaka  (1973). 
For  brevity  we  shall  refer  to  this  model  as  an  "Ohnaka  ramp". 

It  is  used  as  an  example  of  a dislocation  function  in  which 
the  particle  velocity  on  the  fault  is  continuous  in  time. 

On  Fig.  1,  this  model  would  look  identical  to  the  propagating 
ramp.  Fig.  2 shows  the  dislocation  function  3n  time  for  three 
models,  including  the  Ohnaka  ramp  which  was  omltu  u from  Fig.  1. 

For  all  the  models  given  below,  rupture  begins  at  a point 
or  on  a line  at  t = 0,  and  a given  point  on  the  fault  plane 
ruptures  at  t = tr(£j,  £2)*  The  dislocation  direction  is 
always  parallel  to  Ci:  D * (D,  0,  0). 

The  point  t - (0,  0)  is  taken  to  be  the  geometrical  center 
of  the  fault  for  all  models  except  the  pinned  annealing  crack, 
where  it  is  the  point  rupture  initiates. 
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Propagating  Ramp 
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Elliptic  Propagating  Ramp 
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Self  Similar  Crack 
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Bilateral  Propagating  Ramp 
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Numerical  Results 


To  permit  meaningful  comparisons  between  different  models, 
two  normalizations  are  used  through  the  rest  of  this  paper 
unless  specifically  noted  otherwise.  First,  „e  take  the  fault 
area  km2 . This  fault  area  might  be  in  the  range  appro- 

priate for  a magnitude  2',  to  3 earthquake  (Wyss  and  Brune, 
1968).  Second,  we  note  from  equation  (1)  and  from  liaskell 
(1969)  that  the  displacement  at  the  station  is  linearly  pro- 
portional to  the  displacement  on  the  fault  surface ..  Therefore , 
without  loss  of  generality  we  take  the  average  displacement 

fault  to  be  u - 1 cm.  The  displacements  calculated 
are  plotted  as  the  fraction  of  the  average  displacement  on 

the  fault.  Unless  noted  otherwise,  the  station  is  a distance 
1 km  from  the  center  of  the  fault,  but  may  be  considerably 
closer  to  some  portion  of  the  fault. 

In  what  follows,  all  the  above  dislocation  functions  will 
be  Uoed  to  calculate  displacements  at  several  stations.  When 
comparing  faults  with  different  geometry,  the  parameters 
are  chosen  and  the  final  faults  aligned  so  the  models 

have  a.  much  of  the  fault  plane  in  common  as  possible. 

The  location  of  all  the  stations  used  are  listed  in  Table  i. 

We  will  be  examining  the  derived  displacements,  and  comparing 
to  see  What  can  be  learned  about  the  source  model.  The  crucial 
Point  in  examining  the  differences  between  two  seismograms 
see  if  tho„e  differences  are  large  enough  to  make 


it  possible  to  distinguish  between  the  two  source  models. 

In  many  cases,  we  find  two  different  source  models  are  indis- 
tinguishable, iri  that  their  near  field  motions  differ  by  less 
than  the  uncertainty  found  in  comparable  real  data. 

Fig.  3 shows  the  results  of  calculations  for  three  differ- 
ent models  of  propagating  ramps.  The  three  models  are 

1.  a starting  model  of  v = 3.0  km/sec,  t = 0.12  sec, 

2.  a model  where  t has  been  decreased  to  0.0*155  sec, 

3.  a model  with  the  lower  value  of  r , and  where  the 

rupture  velocity  has  been  decreased  to  2.5  kn/sec. 

It  is  seen  in  Fig.  3a  that  models  (1)  and  (3)  look  remarkably 

similar.  Fig.  3b  shows  results  of  the  same  models  at  another 

station,  with  the  same  results.  These  figures  illustrate  the 

important  point  that,  when  a propagating  ramp  is  chosen  to 

fit  a displacement  record,  the  choices  of  rupture  velocity 

and  rise  time  can  be  traded  off. 

It  is  interesting  for  these  propagating  ramp  models  to 

compare  the  total  time  required  for  rupture.  V.'e  designate 
as  Tr  the  total  time  interval  during  which  the  fault  is 

rupturing.  For  a propagating  ramp,  Tr  = L/v  + t.  For  the 

cases  shown  in  FLg.  3a  and  3b, 

1.  T ••=  + °-12  “ -lt43  sec. 

r 3 - u 

2 T = + 0.0*155  ” .3^8  sec. 

r 3.0 

3.  Tr  = -°-^p  + 0.0*155  = .^33  sec. 

cases  which  resemble  each  other  closely,  Tr  differs 


For  the 


by  about  2*.  For  the  cases  which  appear  significantly  differ- 
ent, Tr  differs  by  about  205.  This  suggests  that  Tr  is  an 

important  parameter  in  controlling  strong  motion  near  a fault. 
To  investigate  this  further,  additional  cases  were  run.  For 
each  of  these  cases,  two  propagating  ramp  models  were  chosen 
such  that  the  rupture  velocities  differed  substantially  but 

Tr  was  the  same.  In  Fig.  3c.  the  displacements  are  shown  at 
seven  stations  representing  several  azimuths  from  the  center 
of  the  fault.  Thus  we  are  testing  whether  the  agreement  m 
Pies.  3a  and  3b  Is  a result  of  the  choice  of  stations.  It  can 
be  seen  that  the  waveform  varies  strongly  with  azimuth  where 
bhe  rupture  Is  propagating  toward  the  station,  but  fairly 
weakly  with  azimuth  at  back  angles.  However,  at  each  station, 
the  two  propagating  ramps  generate  similar  waveforms.  This 
eire,  therefore.  Indicates  that  the  ambiguity  of  rupture 
velocity  and  rise  time  observed  In  Figs.  3a  and  3b  was  not 
peculiar  to  the  choice  of  station. 

For  Fig.  3c-.,  the  same  rupture  model5  ^ ^ ^ ^ pig_ 

3c,  but  the  station  is  U km 

q k fl0ni  the  center  of  the  fault.  The 

atation  is  near  nodal  for  P-waves.  The  strong  S-„ave  pulses 
two  models  are  quite  similar.  This  example  suggests 
that  within  the  near  field  the  ambiguity  is  not  sensitive  to 
the  distance  between  the  fault  and  the  station. 

Pop  Fig.  3e,  rather  than  a nearly  square  rupture,  the 

fault  has  L - 2.5  km  and  W . CM,  km.  Again,  with  rupture 

toward  the  station  and  m 

' and  Tr  the  same  for  the  two  models,  the 

displacements  derived  from 

from  the  two  models  are  remarkably  similar, 


-i.n  Pigs.  ^-6,  v.'e  show  the  results  of  computations  to  exa- 
mine other  aspects  of  rupture  for  models  which  are  basically 
unilateral.  The  parameter  Tr  is  similar  for  models  being 
compared  in  nearly  all  following  figures. 

Figs,  and  ^b  show  1 1 10  comparison  of  a propagating  ramp 
and  an  Onnaka  ramp.  At  the  two  stations  considered,  the  dis- 
pla  .ement-j  resulting  from  the  two  dislocation  functions  are 
almost  indistinguishable . 

Next  we  consider  tnree  rupture  models  with  slightly  differ- 
ing  fault  ouit aces . Fig.  5 shows  the  faults  used  for  a propa- 
gating ramp , elliptical  propagating  ramp,  and  a pinned  annealing 
crack.  Fig.  o shows  the  resulting  displacement  seismograms 
for  one  station  near  the  initiation  of  rupture  (Pig.  6c)  and 

two  near  where  it  finally  ceases  (Figs.  6a,  6b).  For  this  and 
later  figures  where  the  epicenters  of  the  models  are  not  the 
same,  we  have  shifted  some  of  the  mocels  slightly  along  the 
time  axis.  The  effect  of  t.nis  shift  is  to  have  the  rupture 
(or  in  one  case  later  the  stopping)  occurring  at  about  the  same 
time  and  place  for  all  models.  Apparently  the  effect  of  the 
shape  of  the  fault  is  not  of  major  importance,  since  the  two 
ramps,  with  Identical  dislocation  functions,  differ  only 
slightly.  The  models  show  the  greatest  differences  near 
where  the  rupture  ceases. 

The  difference  between  the  propagating  ramp  and  the  pinned 
annealing  crack  appears  to  be  of  the  nature  that  can  be  com- 
pensated by  a decreased  rupture  velocity  or  increased  rise 


time  in  the  propagating  ramp,  as  in  PigS.  3a>  3b.  However> 
continued  searches  showed  that  the  propasating  ramp  parameters 
used  in  Pig.  6 cannot  be  substantially  Improved  upon,  within 
the  constraint  that  the  average  displacement  be  1 cm.  The 
differences  appear  to  arise  from  the  very  high  particle 
velocities  on  the  fault  surface  of  the  pinned  annealing 
crack.  Since  the  amplitude  of  far-field  (1/r)  radiation 

from  a fault  element  is  proportional  to  the  particle  velocity 

on  that  fault  elppprii-  *.u  . 

ament,  there  is  very  high  radiation  from  the 

crack  tip.  In  fitting  with  a propagating  ramp  model,  this 

n.ay  result  in  as  si  mi  nr  n 

fgnlng  a larger  average  displacement  to  the 

fault  surface  than  actually  occurred. 

So  far,  then,  we  have  examined  four  models  of  rupture 
which  are  predominantly  unilateral  in  nature,  it  appears 
that  at  least  at  these  distances,  and  for  the  fault  sites 

, the  pa.  tide  notion  would  be  very  difficult  to  learn 
from  displacement  records  of  the  earthquake.  The  addition 

of  new  stations  at  the  same  distance  provides  little  new 
information. 

We  now  turn  our  attention  to  two  basically  bilateral  source 
unctions  and  begin  by  comparing  a biiateral  propagating  ramp 

a unilateral  propagating  ramp  having  the  same  geometry, 
t the  station  shown  in  Pig.  7 there  are  5ubstantlal 

ences  in  the  u3  component  which  would  be  sufficient  to  dis- 

tinguisn  between  the  two  model..  In  Fir:  fl  h 

in  ri6.  8,  however,  a bilatera 

model  with  v = ? o -i  . 

3.0  is  compared  with  a unilateral  model  with 


v =.  6.C.  The  rise  times  are  the  same  Thus  Tp  is  identical 
for  the  two  models.  The  differences  at  the  station  chosen 
in  Pig.  8a  are  small.  Where  the  unilateral  model  is  not  pro- 
pagating predominantly  toward  the  station,  the  two  models 
differ  substantially,  as  shown  in  Pigs.  8b,  8c.  Thus  addition 
of  new  stations  helps  distinguish  between  unilateral  and 
bilateral  rupture.  This  suggests  that  the  predominant  direction 
or  directions  of  rupture  propagation  is  important  i determining 

strong  motion  in  the  near  field. 

_.  t 

Fig.  9 compares  a self  similar  crack  model,  which  grows 

radially,  with  a bilateral  propagating  ramp.  The  bilateral 
crack  grows  with  the  same  propagation  velocity  (in  the 

direction)  as  the  self  similar  crack.  These  cracks  grow  and 
stop  in  very  different  ways,  but  the  predominant  direction 

of  rupture  growth  is  similar  to  the  extent  that  rupture  initiate 
on  a point  or  a line  through  the  center  of  the  fault.  Fig.  9 
shows  differences  in  the  resulting  near-field  displacements. 
However,  these  differences  are  perhaps  not  as  great  as  one 
would  have  expected,  for  such  different  fault  models. 

Although  most  models  of  a unilateral  ramp  show  little  re- 
semblance to  a self-similar  crack,  one  interesting  comparison 
can  be  made.  This  is  shown  in  Pig.  10.  A propagating  ramp 
is  used  with  a propagation  velocity  of  100  km/sec,  so  that 
the  entire  fault  starts  and  stops  nearly  instantaneously. 

There  are  substantial  differences  at  the  beginning  of  the 
records,  reflecting  the  different  way  in  which  the  dislocations 
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begin.  But  the  strong  resenblance  of  the  two  events  at  the 
end  of  the  records  reflects  that  they  stopped  in  a nearly  * 
identical  fashion. 

The  previous  figures  have  not  shown  a particularly  encou- 
raging situation  xor  the  goal  of  identifying  details  of  the 
rupture  process  from  the  displacement  near  the  fault.  It  is 
possible  that  this  is  due,  in  part,  to  the  partic:  ar  dimen- 
sions chosen  in  the  normalisation  described  earlier.  All 
the  earthquakes  shown  so  far  would  probably  correspond  to 
magnitude  around  2^  to  3.  Certainly  as  the  station  is  moved 
closer  to  the  fault  plane,  displacement  there  would  more  closely 
resemble  the  displacement  dislocation  function  for  that  part 
of  the  fault.  One  would  also  expect  that,  if  the  station  were 
kept  the  same  absolute  distance  from  the  fault,  then  for  a 
larger  earthquake,  represented  by  a larger  fault  plane,  the 
displacement  at  the  station  wou  .d  more  closely  resemble  the 
fault  dislocation.  We,  therefore,  imagine  a strike  slip  geo- 
metry with  a station  located  on  the  ground  surface  for  earth- 
quakes with  fault  area  12.6  km2,  70.8  km2  and  *402.0  km2. 

For  convenience,  we  will  label  these  by  magnitude  *4,  5,  and 
7,  respectively,  althcn.  these  fault  areas  are  smaller  than 
the  areas  predicted  for  earthquakes  of  such  magnitudes  by 
the  formula  of  Vryss  and  Brune  (1968).  In  increasing  the  fault 
areas,  we  take  the  predominant  increase  to  be  in  the  length. 

As  a result,  the  length  of  the  "magnitude  7"  earthquake 


developed  for  the  May  18,  1940,  Imperial  Valley,  California, 
earthquake,  which  had  a magnitude  7.1  (Richter,  1958).  The 
stations  are  a distance  5 km  from  the  center  of  the  faults, 
thus  slightly  over  one  fault  length  for  the  "magnitude  V 
earthquake.  We  pose  the  question  of  whether  one  could  dis- 
tinguish between  a propagating  ramp  and  an  Ohnaka  ramp.  Results 
of  these  calculations  are  shown  in  Figs.  11a,  11b,  11c.  The  two 
source  time  functions  give  larger  differences  for  all  three  case 
than  they  did  for  the  smaller  fault  in  Figs.  4a,  4b.  Neverthelo 
it  is  not  clear  that  a distinction  could  be  made  in  any  of 
these  cases,  using  data  with  possible  noise  and  unknown 
surface  effects.  The  components  not  parallel  to  the  direction 
of  dislocation  show  possibly  important  differences,  while  the 
component  parallel  to  the  dislocation  docs  not. 


Discussion 

So  far  we  have  seen  *hat  a considerable  variety  of  dislo- 
cation functions  and  parameters  can  be  used  to  give  basically 
similar  wave  forms  at  stations  near  a moving  fault.  Some 
understanding  of  this  similarity  can  be  obtained  by  examina- 
tion of  the  Fourier  amplitude  spectrum  of  one  of  the  cases. 

v 

This  is  shown  in  Fig.  12  for  one  of  the  propagating  ramp 
models  shown  in  Fig.  3b.  It  should  be  noted  that  this  Fourier 
amplitude  spectrum,  is  not  comparable  to  the  spectra  which  are 
most  often  seen  in  connection  with  source  studies.  There  are 
two  important  differences.  Firs* . it  is  a near  field  spectrum 
rather  than  far  field.  It,  therefore,  behaves  as  at  long 

periods  due  to  the  permanent  offset.  Second,  the  P-  and  S- 

waves  are  not  separated,  and  the  spectrum  includes  both. 


In  Fig.  12  the  displacement  spectral  density  decreases  rapidly 
to  25  Hz,  then  begins  to  rise  again.  This  rise  is  due  to  the 
small  errors  involved  in  the  calculation  of  displacement.  The 
frequency  where  it  appears  is  correlated  approximately  with 
the  average  spacing  of  integration  points  on  the  fault  and 
it  can  be  adjusted  through  the  error  limit  specified  in  the 
calculation . 

If  any  corners  exist  in  the  spectrum,  one  might  expect 
them  to  correlate  with  times  characteristic  of  the  source 
function.  The  rupture  front  crosses  the  fault  in  0.^  sec, 
corresponding  to  a frequency  of  2,5  Hz.  The  rise  time  is 
0.0^55  sec  corresponding  to  a frequency  of  about  22  Hz.  If 
one  assiduously  scare nes  the  spectrum  for  corners,  they  might 
be  found  indicates  at  3*5  Hz  and  12  Hz,  with  a fall-off  of 
w"1*  after  the  second  one.  The  important  point  is  that  before 
the  emergence  of  numerical  noise,  the  spectral  amplitude  is 
decreasing  very  rapidly  with  increasing  frequency.  The  high 
frequencies  which  are  needed  to  examine  details  of  the  fault- 
ing, such  as  how  the  rupture  starts,  were  not  reached  in  this 
calculation  with  assigned  error  of  one-tenth  normal,  or  one- 
half  percent.  These  frequencies  would  be  very  weakly  repre- 
sented in  a record  of  displacement,  even  when  calculated  vjith 
high  accuracy. 

Clearly  it  would  be  advantageous  to  study  the  acceleration 
instead  of  the  displacement.  Since  this  has  the  effect  of 
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multiplying  the  amplitude  spectrum  by  w2  the  acceleration 
response  v;ould  approximately  be  flat,  to  about  20  liz . More- 
over, near-field  data  is  usually  gathered  In  terms  of  accelera- 
tion. Unfortunately,  to  achieve  the  needed  accuracy  for  two 
numerical  aif ferentiations  (as  was  done  in  the  test  for  obeying 
the  wave  equation)  one  requires  at  least  a ten-fold  increase 
in  computation  time. 

It  is  not  at  all  clear  whether  calculations  of  accelera- 
tions are  justified.  Effects  due  to  the  free  surface  and 
due  to  layering  can  be  expected  to  be  large  at  higher  fre- 
quencies. Therefore,  a dislocation  model  formulated  for  an 
infinite  space  would  be  of  questionable  value  in  any  compari- 
son with  data  at  other  than  long  periods.  This  suggests  that 
a dislocation  model,  as  used  here,  is  not  appropriate  for 
learning  about  the  details  of  rupture.  This  does  not  modify 
its  usefulness  in  studying  the  gross  features  of  rupture,  as 
has  been  done  for  the  Parkfield  and  San  Fernando  earthquakes 
in  the  studies  mentioned  earlier. 

An  examination  of  the  displacements  in  Figs.  3-H  will 
show  that  in  many  cases,  for  large  and  small  earthquakes , 
the  component  of  motion  parallel  to  displacement  on  the  fault 
(Ul)  is  the  least  useful  in  distinguishing  between  different 
dislocations.  The  component  perpendicular  to  the  fault  plane 
tends  to  be  the  most  sensitive  to  choice  of  model. 
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Conclusions 


There  Is  considerable  ambiguity  in  using  near-field  dis- 
placements to  determine  parameters  of  a dislocation  model 
for  an  earthquake.  One  fault  length  from  the  center  of  small 
faults,  one  can  determine  quite  well  the  total  tine  spent  in 
rupture.  V/ith  more  than  one  station,  one  can  also  determine 
i-he  predominant  direction(s)  of  rupture  growth.  However, 
with  the  exception  of  very  high  particle  velocities  on  the 
fault  surface,  details  of  the  time  function  or  the  ^ault  shape 
have  relatively  little  influence  on  the  displacements  at  near- 


field stations. 
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Table  1 


Station  Locations 


x j (km) 

x2 (km)  . 

x j (km) 

0.863 

-0.500 

0.075 

0.750 

-0.500 

0.433 

0.433 

-0.500 

0.750 

0.0 

-0.500 

0 . 866 

-O.433 

-0.500 

0.750 

-0.750 

-0.500 

• 

* 0.433 

-O.863 

-0.500 

0.075 

0.837 

-0.500 

0.224 

-0.837 

-0.500 

0.224 

3 • 95*1 

-0.500 

0.346 

0.0 

4.0 

3.0 
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Fig . 1. 


Fig.  2 


Figure  Captions 

Growth  of  crack  models  in  tine.  Solid  lin^s  out- 
line the  area  ruptured  at  each  of  the  tines  shown . 

Broken  lines  are  healing  fronts,  separating  areas 
which  have  ceased  notion  from  areas  which  are  still 
moving.  Names  by  which  these  models  are  referred 
to  are:  1)  propagating  ramp,  2)  elliptical  pro- 

pagating ramp,  3)  self  sinilap  crack,  *1)  pinned 
annealing  crack.  5)  bilateral  propagating  ramp. 

For  the  ramp  models  the  shaded  are^s  are 

in  motion.  For  the  other  two  models,  a shaded 
point  has  achjeved  less  than  one-half  the  final 
average  displacement  of  the  fault  regardless  of 
whether  it  is  stopped  or  moving.  The  symbol  above 
each  model  will  be  used  to  identify  the  model  in 
all  later  figures. 

Particle  motion  on  the  fault  with  time.  To  the 
left,  for  each  of  three  models,  is  an  outline  of 
the  final  fault  shape.  To  the  right  is  the  time 
history  of  particle  motion  for  ^ticles  on  the 
fault,  at  positions  indicated  by  the  beginning 
level  of  the  liner.  The  central  model,  indicated 
by  a new  symbol  in  this  figure,  is  a propagating 
ramp  with  an  Ohnaka  time  function  (see  text). 

For  the  pinned  annealing  crack  at  right  only, 
the  time  functions  apply  only  to  the  midline  of 
the  fault.  The  horizontal  (time)  and  vertical 
(displacement)  scales  are  arbitrary. 
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Effect  of  varying  the  rupture  velocity  and  rise 
time  for  a propagating  ramp.  The  rupture  velocity 
and  rise  tine  for  each  model  is  given  in  the  upper 
right  corner.  Parameters  common  to  all  models 
are  L = 0.968  km,  W = 1.032  km,  Dmax  = 1.0  cm. 

In  the  box  at  lov;er  right,  two  views  of  the  fault 
are  given,  showing  the  fault  and  station  projected 
onto  the  x,  - x3  plane  and  the  x,  - x2  plane 

the  station  is  no.  2 in  Table  1.  For  the  view 
in  the  x 1 - x 3 plane,  looking  in  the  +x2  direction. 
The  arrows  indicate  that  the  offset  is  left  lateral. 
Fault  growth  is  in  the  same  sense  as  in  Fig.  1, 
i.e.,  here  rupture  initiates  at  the  left  and  crosses 
the  fault  to  the  right.  The  two  vertical  lines 
next  to  the  fault  views  show  the  distance,  on 
the  same  scale,  which  would  be  traveled  by  a P- 
wave  (a  = 6.0  km/sec)  and  an  S-wave  (3  = 3.^  km/sec) 
in  one  time  unit.  The  lines  labelled  u,,  u2,  u3 
give  the  displacements  in  the  x,,  x2 , x3  directions 
at  the  station  shown  in  the  box.  The  vertical  scale 
is  the  fraction  of  the  average  displacement  on 
the  fault  plane. 


Fig.  3b. 


Model  same  as  Fig.  3a,  at  station  location  Mo.  3. 


Fie.  3c. 


Comparison  at  seven  stations  of  displacements 
from  two  propagating  ramp  models  with  the  same 


Tv . For  both  models,  L = 1.0  km,  W = 1.0  km, 

Dn,ax  ~ 1.0  cm.  The  rupture  velocities  and  rise 
times  are  indicated.  The  numbers  beside  the 
stations  refer  to  the  station  number  in  Table  1. 

Fig.  3d.  Comparison  of  the  models  in  Fig.  3c  at  station 
10,  lJ.0  km  from  the  center  of  the  fault. 

» 

Fig.  3e . Comparison  of  two  long,  narrow  propagating  ramps 
with  identical  T . Model  parameters  are  L = 2.5 
km,  V/  = O.lJ  km,  Dmax  = 1.0  cm.  Rupture  velocity 
(in  km/sec)  and  rise  time  (in  sec)  for  each  model 
is  shown  in  the  figure  Station  No.  8.  For  explana- 
tion of  other  symbols,  see  Fig  3a. 

Fig.  JJa,b.  Comparison  of  propagating  ramps  with  a linear 
and  a smoothed  time  function  at  Station  No.  2 
(*Ja)  and  3 (^b).  Model  parameters  are  L = O.968 
km,  W = 1.032  km,  v = 3.0  km/sec,  Dmax  = 1.0  cm 
for  both  models.  For  the  linear  time  function, 
x = 0.12  sec,  and  for  the  Ohnaka  time  function 
x q = 0.03  sec.  For  explanation  of  other  symbols, 
see  Fig.  3a. 


Fig.  5. 


Comparison  of  the  fault  rupture  areas  of  the 
models  in  Figs.  6.  The  three  stations  used 
are  shown  projected  onto  the  fault  plane.  Arrows 
point  to  where  rupture  initiates. 


Comparison  of  three  faulting  models  at  three 
stations.  Parameters  for  the  models  are  the 
following:  Propagating  ramp:  L = 2.31  km, 

V/  = 0.^3  km,  v = 3.0  km/sec,  t " 0.0^5  sec, 

* 

Dmax  = 1.0  cm.  Elliptical  propagating  ramp: 

L = 2 . 60  km,  W = 0.^9  km,  v = 3.0  km/sec, 
t = 0.05  sec,  Dmax  = i.o  cm.  Pinned  annealing 
crack:  s = 6.17  cm/sec,  a = 2.7  km/sec, 

v = 1.0  km/sec,  a}  = -0.191  kn,  a2  = 0.0  km, 
c = 3.2  km/sec.  Stations  are  No.  8 (6a), 

No.  2 (6b),  and  No.  9 (6c).  For  explanation 
of  other  symbols,  see  Fig.  3a. 


Comparison  of  a unilateral  and  a bilateral 
faulting  model.  For  loth  models,  L = 1.0  km 
VI  = 1.0  km,  v = 2.5  km/sec,  t = 0.05  sec, 

Dmax  = 1*0  cn  Station  No.  2.  For  explanation 
of  other  symbols,  see  Fig.  3a. 
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Pig.  8a,b,c.  Comparison  of  a unilateral  and  a bilateral 

faulting  model  at  tv/o  stations.  For  both 
models,  L = 1.0  km,  W = 1.0  km,  t = 0.2  sec, 

Dmax  = 1*0  cm.  For  the  bilateral  model,  v = 

3.0  km/sec,  and  for  the  unilateral  model, 
v = 6.0  kn/sec.  Stations  are  I Jo . 2 (3a), 

No.  9 (8b),  and  No.  ^ (8c).  For  explanation 
of  other  symbols,  see  Fig.  3a. 

Fig.  9a,b,c.  Comparison  of  a bilateral  and  a sel^f  similar 

crack  model  at  three  stations.  For  the  Di- 
lateral crack,  L = 1.0  km,  W = 1.0  km,  v = 

3.0  km/sec,  t = 0.05  sec,  =1.0  cm:  for 

the  self  similar  crack,  s = ^.125  cn/sec, 
o = 3-0  km/sec,  v = 3.2  km/se<- , t „ v = 0.132 
sec.  Stations  are  No.  2 (9a),  No.  3 (9b), 
anc  4 (9c).  For  explanation  of  other 
symbols,  see  Fig.  3a. 

Fig.  10a, b,c.  Comparison  of  a self  similar  crack  and  a pro- 
pagating ramp  which  stops  similarly.  For  the 
self  similar  crack,  s = A. 125  cm/sec,  a = 

3.0  km/sec,  v = 3.2  km/sec,  t = 0.182  sec; 

rndx 

for  the  propagating  ramp,  L = 0.968  km,  VJ  = 

1.032  km,  v = 100.0  km/sec,  t = 0.182  sec, 

Dnax  = 1.0  cm.  Stations  are  No.  2 (10a),  No. 

3 (30b),  and  No.  4 (10c).  For  explanation  of 
other  symbols,  see  Fig.  3a. 


Fig.  11a. 


Comparison  of  propagating  ramps  with  linear 
and  smoothed  time  functions.  The  parameters 
might  be  appropriate  to  an  earthquake  of  mag- 
nitude about  4.  For  both  models,  L = 3.55 
km,  W = 3-55  km,  v = 3.0  km,  Dmax  = 1.0  cm. 
For  the  linear  ramp,  x = 2.5  sec.  For  the 
smoothed  ramp,  tq  = 0.625  sec.  In  this  case, 
the  station  No.  11  is  a distance  5.0  km  from 

the  center  of  the  fault.  For  explanation  of 
other  symbols,  see  Fig.  3a.  — 

Fig.  lib.  A "magnitude  5"  earthquake.  For  both  models, 

L = 11.8  km,  W = 6.0  km,  v = 3.0  km/sec, 

Dmax  = 1.0  cm.  For  the  linear  ramp,  t = 5.0 
sec.  For  the  smoothed  ramp,  t0  = 1.25  sec. 
Station  No.  11. 

Fig.  11c.  A "magnitude  7"  earthquake.  For  both  models 

L = 6 7.0  km,  W = 6.0  km,  v = 3*0  km/sec, 

Dmax  = 1.0  cm.  For  the  linear  ramp,  t = 5.0 
sec.  For  the  smoothed  ramp,  t0  = 1.25  sec. 


Station  No.  11. 


Fig.  12. 
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Fourier  Displacement  spectral  amplitude  for 
one  of  the  models  shown  in  Fig.  3b.  This  is 
the  spectrum  of  near  field  displacement,  so 
the  long  period  limit  should  be  as  of1.  Note 
that  the  spectrum  contajns  both  P-  and  S-waves. 


AMPLITUDE  SPECTRUM 


